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We study the superconducting state in recently discovered high-T c superconductor Ka;Fe2Se2 based 
on the ten-orbital Hubbard-Holstein model without hole-pockets. When the Coulomb interaction is 
large, spin- fluctuation mediated d-wave state appears due to the nesting between electron-pockets. 
Interestingly, the symmetry of the body-centered tetragonal structure in K x Fe2Se2 requires the 
existence of nodes in the d-wave gap, although fully-gapped d-wave state is realized in the case of 
simple tetragonal structure. In the presence of moderate electron-phonon interaction due to Fe- 
ion optical modes, on the other hand, orbital fluctuations give rise to the fully-gapped s++-wave 
state without sign reversal. Therefore, both superconducting states are distinguishable by careful 
measurements of the gap structure or the impurity effect on T c . 

PACS numbers: 74.20.-z, 74.20.Fg, 74.20.Rp 



The pairing mechanism of high-T c iron-based super- 
conductors has been significant open problem. The 
main characters of FeAs compounds would be (i) the 
nesting between electron-pockets (e-pockets) and hole- 
pockets (h-pockets), and (ii) the existence of orbital de- 
gree of freedom. By focusing on the intra-orbital nest- 
ing, fully-gapped sign-reversing s-wave state (s±-wave 
state) had been predicted based on the spin fluctuation 
theories^* 2 -. On the other hand, existence of moderate 
electron-phonon (e-ph) interactions due to Fe-ion opti- 
cal phonons and the inter-orbital nesting can produce 
large orbital fluctuations 3 -. Then, orbital-fluctuation- 
mediated s-wave state without sign reversal (s ++ -wave 
state) had been predicted by using the random-phase- 
approximation (RPA) 3,4 or the fluctuation-exchange 
(FLEX) approximation 5 . According to the analysis in 
Refs. [y0 the s ++ -wave state is consistent with the 
robustness of T c against randomness^ as well as the 
"resonance-like" hump structure in the neutron inelas- 
tic scattering 1 —. Non-Fermi liquid transport phenomena 
in pii can be explained by the development of orbital 
fluctuations 5 -. 

Recently, iron-selenium 122-structure compound 
A x Fe2Se2 (A= alkaline metals) with T c ~ 30 K was 
discovered 12 . This heavily electron-doped superconduc- 
tor has been attracting great attention since both the 
band calculation s 13 ! 14 and angle-resolved-photoemission- 
spectrum (ARPES) measurements^ 5 - - — indicate the 
absence of h-pockets. NMR measurements reports the 
weakness of spin fluctuations^, and both ARPES^ 5 - - — 
and specific heat measurements 19 indicate the isotropic 
SC gap. Thus, study of A ;E Fe2Se2 will give us important 
information to reveal the pairing mechanism of iron 
pnictides. 

The unit-cell of iron-based superconductors contains 
two Fe atoms. However, except for 122-systems, one 
can construct a simple "single-Fe model" from the origi- 
nal "two-Fe model" by applying the gauge transforma- 
tion on d-orbitals 20 . By this procedure, the original 



Brillouin zone (BZ) is enlarged to the "unfolded BZ". 
Based on the single-Fe model, spin-fluctuation-mediated 
d-wave state (B\ g representation) "without nodes" had 
been proposed 2 ^"—, by focusing on the nesting between 
e-pockets. However, we cannot construct a "single-Fe 
model" for 122 systems since finite hybridization between 
e-pockets prevents the unfolding procedure 2 ^. Therefore, 
theoretical study based on the original two-Fe model is 
highly desired to conclude the gap structure. 

In this paper, we study the ten-orbital (two Fe atoms) 
Hubbard-Holstein (HH) model for KFe2Se2 using the 
RPA. When the Coulomb interaction is large, we ob- 
tain the (i-wave SC state due to the spin fluctuations, as 
predicted by the recent theoretical studies in the single- 
Fe Hubbard models 2 ^ - —. However, the gap function on 
the Fermi surfaces (FSs) inevitably has "nodal structure" 
in the two-Fe model, due to the symmetry requirement 
of the body-centered tetragonal lattice. On the other 
hand, orbital-fluctuation-mediated s ++ -wave state is re- 
alized by small e-ph coupling; the dimensionless coupling 
constant A = gN(0) is just ~ 0.2. Since the nodal SC 
state is fragile against randomness, study of impurity ef- 
fect will be useful to distinguish these SC states. 

We perform the local-density-approximation (LDA) 
band calculation for KFe2Se2 using Wien2k code based 
on the experimental crystal structure^ 2 -. Next, we de- 
rive the ten-orbital tight-binding model that reproduces 
the LDA band structure and its orbital character using 
Wannier90 code and Wien2Wannier interface 24 . The dis- 
persion of the model and the primitive BZ are shown 
in Figs. [1] (a) and (b). Based on a similar ten-orbital 
model, Suzuki et al studied the s±-wave gap structure 
for BaFe 2 As 2 — . 

In Fig. m we show the FSs of KFe 2 Se 2 for (c) k z = 
and (d) k z = tt planes when the electron number per Fe- 
ion is n = 6.5: On each plane, there are four large and 
heavy e-pockets around X and Y points, and one small 
and light e-pockets around Z point. For n = 6.5, the 
energy of the h-band at T point from the Fermi level, 
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FIG. 1: (color online) (a) Dispersion of the present ten-orbital 
model for KFe2Se2- T, X, and Z points are on the k z = 
plane, (b) Primitive BZ for body-centered tetragonal lattice. 

(c) (d) FSs on the k z = plane and k z — tt plane. The green 
(light gray), red (gray), and blue (dark gray) lines correspond 
to xz, yz, and xy orbitals, respectively. The diamond-shaped 
shadows in (d) indicates the sign of basis function for B\ g 
representation, (e) Weight of each orbital on the inter FS 
(FS1) and the outer FS (FS2) as function of 0; is shown in 

(d) . (f) Hybridization between two e-pockets in 122 systems 
should create the nodal d-wave gap. 



Eh, is about —0.07 eV. Since the obtained FS topology 
and the value of Eh are consistent with recent reports by 
ARPES measurements^ - —, we study the case n = 6.5 
hereafter. In the present BZ in (b), V and Z points and 
X and Y points in (c) are not equivalent, and k z = 2tt 
plane is given by shifting (c) by (7r,7r). As for (d), T and 
T' points and P and P' points are equivalent, meaning 
that the reciprocal wave vector on the k z = tt plane is 
(7r,7r) and (tt, —tt). The diamond-shaped shadows in the 
k z = tt plane indicates the sign of basis function for B\ g 



(x 2 — y 2 -type) representation, which has nodes on the 
P-P' line on both FS1 (inner FS) and FS2 (outer FS). 

To confirm the existence of nodes, we verify that FS1 
and FS2 in KFe2Se2 are largely hybridized. In fact, the 
weights of d-orbitals on FS1,2 given in Fig. [1] (e) are 
smooth functions of which is the strong evidence for 
the hybridization in wide momentum space. This hy- 
bridization disappears when inter-layer hoppings are ne- 
glected: Then, both xy(F$l) and xy(FS2) show cusps 
at = 7r/4, and xz(FS2) suddenly drops to almost zero 
for 6 > 7r/4. In (f), we explain the origin of nodal gap 
based on the fully-gapped d-wave solution in the single- 
Fe model— - —: By introducing inter-layer hoppings, two 
elliptical e-pockets with positive and negative A in the 
unfolded BZ are hybridized to form FS1 and 2 with four- 
fold symmetry. As a result, nodal lines inevitably emerge 
on FS1 and 2, at least near the \k z \ = tt plane. 

Here, we study the ten-orbital HH model using the 
RPA. As for the Coulomb interaction, we consider the 
intra-orbital term J7, the inter-orbital term U' ', Hund's 
coupling or pair hopping J, and assume the relation U = 
U' + 2 J and J = U/6. In addition, we consider the e-ph 
interaction due to Fe-ion optical phonons; the phonon- 
mediated e-e interaction (—g) and its matrix elements 
are presented in Ref. |4|. Hereafter, we perform the RPA 
on the two-dimensional planes for k z =0, 7r/2, and tt. 

For n = 6.5 and k z = 0, the critical value of g for the 
orbital-density-wave (ODW) is g c = 0.23 eV for U = 0, 
and the critical value of U for the spin- density- wave 
(SDW) is U = 1.18 eV for g = 0. These values change 
only rsj 2% for different k z . The obtained U-g phase di- 
agram is very similar to Fig. 2 in Ref. 4, irrespective of 
the absence of h-pockets in KFe2Se2- The reason would 
be (i) the density-of-states (DOS) in KFe2Se2 is about 
leV -1 per Fe, which is comparable with other iron pnic- 
tides, and (ii) the nesting between e-pockets is rather 
strong because of their square-like shape. Figure [2] (a) 
shows the total spin susuceptibility x s 0) at E7 = 1.1 
eV and g = for k z = plane. x s ls given by the intra- 
orbital nesting, and its peak position is q « (it, Air), 
consistency with previous studies 21 23 . The obtained in- 
commensulate spin correlation is the origin of the d-wave 
SC gap. Figure [2] (b) shows the off-diagonal orbital sus- 
ceptibility Xyztxyiyz&yiQi 0) ^ or the k z = plane at U = 
and g = 0.22 eV; its definition is given in Refs. 3-5. It is 
derived from the inter-orbital nesting between xz and xy, 
and its peak position is q « (0.77T, OAtt). Note that the 
peak position of Xxz,xy,xz,xy is Q ~ (0.4tt, 0.7tt). The ob- 
tained strong spin- and orbital-correlations are the origin 
of the d-wave and s ++ -wave SC states. 

In the following, we solve the linearized gap equation 
to obtain the gap function, by applying the Lanczos al- 
gorithm to achieve reliable results. In the actual cal- 
culation results shown below, we take 64 x 64 fc-point 
meshes and 512 Matsubara frequencies. First, we study 
the spin- fluctuation- mediated SC state for U < U c by 
putting g = 0. Figures [3] (a)- (c) show the gap functions 
of the d-wave solution at T = 0.03 eV for k z = 0, 7r/2, 
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FIG. 2: (color online) (a) X s (q, 0) for U = 1.1 eV and = 0, 
and (b) X^^^Oz, 0) for £7 = and p = 0.22 eV. 



and 7r, respectively. In case of t/ = 1.1 eV, the eigen- 
value Ae is 0.61 for (a), 0.63 for (b), and 0.62 for (c); 
the relation Ae > 1 corresponds to the SC state. They 
are relatively small since the SC condensation energy be- 
comes small when the SC gap has complicated nodal line 
structure. On the (c) k z = tt plane, the nodal lines are 
along = 7r/4 and 37r/4 directions, consistently with the 
basis of Big representation in Fig. [1] (d). These nodes 
move to near the BZ boundary, = and 7r, on the (b) 
k z — k/2 plane, and they deviate from the FSs on the 
(a) k z = plane. As results, the nodal gap appears for 
tt/2 < \k z \ < 3tt/2 in the whole BZ \k z \ < 2tt. 

We also obtain the s±-wave state, with the sign re- 
versal of the SC gap between e-pockets and the "hidden 
h-pockets below the Fermi level" given by the valence 
bands 5,6. The obtained solution is shown in Fig. [3] 
(d) for k z = 7r. Interestingly, the obtained eigenvalue is 
Ae = 0.99 for U = 1.1 eV, which is larger than Ae for 
d-wave state in Fig. E9(a)-(c). Such large Ae originates 
from the scattering of Cooper pairs between e-pockets 
and the "hidden h-pockets" , which was discussed as the 
"valence-band Suhl-Kondo (VBSK) effect" in the study 
of Na^CoOs in Ref. HE 

Here, we analyze the T-dependence of Ae based on a 
simple two-band model with inter-band repulsion: The 
set of gap equations is given by 26 AeA^ = —VN e L e A e 
and AeA 6 = —VN^LhAh, where V > is the repul- 
sive interaction between e- and h-pockets, and N Gj h is 
the DOS near the Fermi level. When (i) the top of 
the h-pocket is well above the Fermi level, L e = Lh = 
ln(1.13co> c /T), where uo c is the cutoff energy. Thus, the 
eigenvalue is given as Ae = VyJN e Nu ln(1.13o; c /T) oc 
— InT, similar to single-band BCS superconductors. On 
the other hand, when (ii) h-pocket is slightly below the 
Fermi level, L h = (1/2) \n(uj c /\E h \), where E h < is the 
energy of the top of h-band 26 . Thus, the eigenvalue is 
given as A E = Vy/N e N h L hy /]n(1.13LJ c /T) ex yj- InT. 
Therefore, in case (ii), the T-dependence of Ae is much 
moderate. In fact, as shown in Fig. 31(a), Ae for <i-wave 
state increases monotonically with decreasing T, while 
Ae for s±-wave state saturates at low temperatures. This 
result suggests that the d-wave state overcomes the s±- 
wave state at T c ~ 30K in K x Fe2Se2. Although T c in the 
s±-wave state is ~ 0.06 eV in Fig. HI (a), it is greatly re- 
duced by the self-energy correction that is absent in the 
RPA£. 
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FIG. 3: (color online) SC gap functions for outer FS (FS2) on 
band 7 and inner FS (FS1) on band 8: d-wave gap functions 
on the (a) k z = 0, (b) k z = tt/2, and (c) k z = tv planes. Black 
and green (gray) lines represent the FSs and gap nodes, (d) 
s±-wave gap function on the k z = n plane. 



We discuss the VBSK effect for s± wave state in more 
detail: According to inelastic neutron scattering mea- 
surement of Ba(Fe,Co)2As2 10 , the characteristic spin- 
fluctuation energy is cj s f ~ 100K just above T c ~ 30K. If 
we assume a similar cj s f in KFe2Se2 since T c is close, 
we obtain the relation uj c ~ oj s f <C \Eh\ in KFe2Se2- 
Since Lh is a monotonic decrease function of \Eh\/uj c 
and Lh < 1 for —Eh/u c > 0.15, we consider that d- 
wave state overcomes the s±-wave state in KFe2Se2, as 
far as the spin-fluctution mediated superconductivity is 
considered. Although high-T c s±-wave state mighe be 
realized for \Eh\/uj c < 0.1, then the realized T c will be 
very sensitive to Eh or the filling n£&. 

Now, we study the s ++ -wave state due to orbital fluc- 
tuations on the k z = plane with n = 6.5. In Fig. [H 
(b), we show the a c -dependence of Ae at T = 0.03 for 
the s ++ -wave state with [7 = 0, and the a s -dependence 
of Ae for the d-wave state with g = 0. Here, a c (a s ) 
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phonon energy; ujd =0.15 eV, considering that the calcu- 
lating temperature is about ten times larger than the real 
T c . The SC gap functions for s ++ -wave state are rather 
isotropic, as shown in Fig. |U(c). However, the obtained 
SC gap becomes more anisotropic in case of U > C£. 

We stress that the RPA is insufficient for quantitative 
study of Ae since the self-energy correction E is dropped: 
In Ref. [5|, we have studied the present model based on 
ther FLEX approximation, and found that the critical 
region with a c > 0.95 is enlarged by the inelastic stac- 
cering 7 = ImS. Also, the 7-induced suppression in Ae 
for d- or s±-wave states is more prominent than that for 
5 ++ -wave state, since 7 due to spin fluctuations is larger 
than that due to orbtial nuctuaiotns 5 . 

Recently, we found the paper by Mazin 27 , in which 
the k z dependence of the nodal d-wave gap in Fig. 4 
corresponds to Figs. [3](a)-(c) in the present work. 

In summary, we studied the mechanism of supercon- 
ductivity in KFe2Se2 based on the ten-orbital HH model 
without h-pockets. Similar to iron-pnitide supercon- 
ductors, orbital-fluctuation-mediated s ++ -wave state is 
realized by small dimensionless e-ph coupling constant 
A = gN(0) ~ 0.2. We also studied the spin- fluctuation- 
mediated d-wave state, and confirmed that nodal lines 
appear on the large e-pockets, due to the hybridiza- 
tion between two e-pockets that is inherent in 122 sys- 
tems. Therefore, careful measurements on the SC gap 
anisotropy is useful to distinguish these different pair- 
ing mechanisms. Study of impurity effect on T c is also 
useful since d-wave (and s±-wave) state is fragile against 
impurities. 
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FIG. 4: (color online) (a) T-dependence of Ae for d- and s±- 
wave states. Ae at T = 0.01 eV is underestimated because 
of the shortage of k- and Matsubara-meshes. (b) a s - (« c -) 
dependence of Ae for d-wave (s++-wave) state at T = 0.03 
eV. (c) SC gap functions for s++-wave state. 



is the charge (spin) Stoner factor introduced in Ref. U 
a c = 1 (a s = 1) corresponds to the ODW (SDW) state. 
In calculating the s ++ -wave state, we use rather larger 
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